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^ ' Abstract 

We investigate the adjoint SU(2) lattice gauge theory in 3+1 dimensions with the Wilson pla- 
quette action modified by a Z2 monopole suppression term. For the zero-twist sector we report on 



■ indications for the existence of a finite temperature transition decoupled from the unphysical bulk 
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I. INTRODUCTION 



Pure 3+1 dimensional SU{N) lattice gauge theories in the fundamental representation 
show a finite temperature deconfinement phase transition [1,^ together with the sponta- 

nn 

neous breaking of a global center symmetry p|, |J| governing the critical indices, which 
e.g. for SU(2) correspond to the universality class of the 3D Ising model P]. Lattice uni- 
versality arguments are commonly used to claim that the same should hold for any possible 
lattice action discretization, in particular with different gauge group representations. On 
the other hand, if confinement of quarks can be ascribed to peculiar (topological) exci- 
tations of the continuum Yang-Mills fields, it is not clear how center symmetry breaking 
can lead to an effective theory of QCD At finite temperature vortices classified through 
■Ki{SU{N) /Ijn) ~ along a compactified dimension can of course provide a suitable degree 
of freedom for the pure Yang-Mills case |7|1, but their fate in the presence of fermions re- 

riQ 

mains challenging. Moreover, in SUSY Yang- Mills theory j8|, |9|] and in the Georgi-Glashow 
model confinement is driven by magnetic^ monopoles. Recent work regarding Yang- 



Mills theories based on exceptional groups 



111 ] conjectures that vortices might not be at 



all necessary to have a confining theory. What role do vortices play then in SU{N)1 An 
investigation of the SU (2) lattice gauge theory in the center blind adjoint representation, 
i.e. 5*0(3) ~ SU{2)/Z2, might offer some interesting insight. The presence of finite tem- 
perature effects for such a model has been debated for a long time. Two main problems 
have been faced. First of all Polyakov's center symmetry breaking mechanism is available 
only for the half-integer representations. For integer ones the Z2 local invariance makes the 
fundamental Polyakov loop of no use. The adjoint Polyakov loop, on the other hand, cannot 
be strictly speaking an order parameter for a transition, assuming it existed, since gluons 



will screen adjoint quarks at some distance (~ 1.25 fm for SU{2)) [12, ll3|. However, it 
can still be taken as a signature to distinguish, at finite temperature, a confined although 
screened phase from a deconfined one. Anyhow, the problem of establishing a well defined 
criterion for confinement in this case remains an open and interesting one. 

Second, and most importantly, lattice artifacts lead to first order bulk phase transitions 
at strong coupling, preventing the continuum limit to be reached within the ordinary con- 



ined phase 



I4L Il5| . Similar phase diagrams are shared by SU (N) theories with > 3 



lfi| . For SU{2) 1,2 monopoles were observed to drive the bulk transition 12, Another 
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interesting development came as the vortex free energy, measured on the lattice already in 
the fundamental theory was recently computed also in the adjoint theory [2^, thus 
suggesting how center vortices and twist sectors are entering a center-blind theory. The lat- 
ter work, however, falls short of any attempt to investigate the theory at finite temperature, 
stumbling on ergodicity problems of the algorithm at the bulk phase transition, where tun- 
nelling among different twist sectors becomes strongly suppressed. Thus, a thorough finite 
temperature investigation of the 5*0(3) theory taking into account the topological excita- 
tions is still missing. Several attempts searching for a decoupling of the finite temperature 
transition from the bulk transitions were originally undertaken by S. Datta and R. Gavai 



see 
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23[ and further citations therein). These authors used the fundamental-adjoint 
mixed action modified by Z2 monopole and vortex suppression terms within the Villain-type 
formulation. By studying the specific heat a finite temperature transition consistent with 
the Ising universality class was found. In the pure center-blind adjoint case including only 
monopole suppression it could however only be established for the smallest time-extension 
Nt = 2. Moreover, the relevance of different twist sectors had not yet appeared in the 
literature. 

Here we want to go a step further employing Wilson's adjoint action formulation modified 
by a similar Z2 monopole suppression. In this case the action itself is manifestly center-blind. 
The corresponding phase structure looks in many respects similar to the Villain case, but 
it differs - as we shall show - substantially in the fundamental-adjoint coupling plane. The 
areas along the two axes are completely separated by a bulk phase transition and thus the 
proof of universality - if the latter is really fulfilled - becomes more complicated. In our 
investigation the twist variables will prove an important ingredient to understand the phase 
structure of the model. 

The outline of the paper is as follows. In Section |Tll we introduce the model and discuss 
its phase structure at vanishing temperature with standard bulk observables. In Section 
twist variables will be used in order to show that the first oder bulk transition seems to 
weaken to second order in a certain — A range. In Section HVl we report on the existence 
of the finite temperature transition and check scaling for the critical temperature. Section 

S contains our conclusions. Reports of this work at early stages have been published in 
,|25|. 
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II. ADJOINT ACTION MODEL WITH %2 MONOPOLE SUPPRESSION 



We study the SU{2) lattice gauge theory with a mixed fundamental-adjoint representation 
Wilson action modified by a chemical potential term suppressing Z2 monopoles 

s = E [t^- (1 - ^) + (1 - ^) ] + ^ Ed - ' (1) 

where cTc = Ylpedc^^S^C^^F^p) taken as a product around elementary 3-cubes c defines 
the Z2-charge. For these monopoles a density can be introduced as 

^ = i-<7^E^^) (2) 

normalized such that it tends to unity in the strong coupling region {Nc denotes the number 
of 3-cubes on the 4D lattice). We will be particularly interested in the pure adjoint theory, 
i.e. = 0. In the latter case one can analyze the model with the link variables represented 
both by 50(3) or SU (2) matrices, exploiting the property Tr^ = Tr^ — 1 for the Wilson term 
or picking a random SU{2) representative of the S0{3) link to construct the Z2 monopole 
contribution. As expected nothing changes in the phase diagram, the integration over the 
fundamental links simply doubling the integration domain in the partition function. A 
standard Metropolis algorithm has been used to update the links in both cases and we use 
either one according to the best performance in the case at hand. Although cXc is constructed 
in terms of fundamental representation plaquettes, it is a natural 5*0(3) quantity. In fact, 
for every given 50(3) link variable, the corresponding SU{2) representative can always be 
determined up to a sign. But the latter is cancelled in the product of plaquettes around any 
3-cube, since each link occurs twice for each edge. In other words 

U^{x) ^ -U^{x) ac^ac, V/i, x,c. (3) 

This ensures that for /^^ = the action we study is center-blind in the entire /3a — A plane. 

Previous finite temperature investigations were mainly done in the Villain discretization 
for the 50(3) term in the action, introducing an independent Z2-valued plaquette variable 
23I l2q . Two terms with chemical potentials were added in order to suppress 



CTp 



completely the effect of the lattice artifacts, Z2 monopoles and vortices, whose densities in 
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FIG. 1: The bulk phase transition of the mixed fundamental-adjoint action model ^ with 
monopole suppression for A = 1.0. 

this case are given by 

^=1-<]^E^^)' ^c=llap, (4) 

^ C PGdc 

' P£dl 

where c and / label the Nc 3-cubes and the Ni products of plaquettes having a link in 
common, respectively. We can see that the monopole suppression term in Eq. Q looks 
formally identical with the one used in connection with the Villain-type action, but its 
realization is different and leads to a different phase structure. The Villain discretization 
DC proven to be equivalent to the standard fundamental Wilson action for 
2q (see Eq. ©). On the other hand one can show that the limit A ^ oo 
of the Wilson adjoint formulation can be mapped, in the trivial twist sector, to some positive 
plaquette model . Such a mapping however, as we shall see below, is not equivalent to a 
continuous connection with the standard fundamental action theory. A simple inspection of 
the Pa{v) — Pf phase diagram in these two limits shows that the differences are conspicous. 



fon W^) can 



Ay — > cxo j2^. 
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25l |. One should always bear in 



Therefore, Ay and A must not simply be identified 
mind that in the Villain case the S0{3) invariance under U^{x) — >■ —U^{x) is not realized 
in the action itself but only once the auxiliary crp variables are integrated out. In Fig. ^ 
the (3f — (3a phase diagram of our model (j2I) at T = is shown for rather strong monopole 
suppression A = 1.0. Phase I is connected with the ordinary confinement phase of the 
standard SU (2) model within the fundamental representation, whereas phase II extending 



to /3a oo occurs completely decoupled from phase I for > 0. Indeed, simulating 
the model with the lattice size up to 12^ the average plaquette variable exhibits a strong 
discontinuity across and metastable states on top of the transition line. At (3a = 0.0 the 
latter was located at (3f — 0.96. By studying the fundamental representation Polyakov loop 
and its susceptibility for lattice size 4 x 12^ the finite-temperature phase transition was seen 
at Pf — 1.35, i.e. within phase I as one would expect. The phase structure clearly differs 
from that of the Villain case, where the transition line for increasing (3f extends along the 
latter axis up to — > oo (see jbsl). Therefore, universality of physics observed for (3f = ^ 
in comparison with (3a = ^ remains an open question. 

3 
2.5 
2 

^ 1.5 
1 

0.5 


0.2 0.4 0.6 0.8 1 1.2 1.4 







1^' order bulk i — 
2"" order bulk ^-X-^ 






n 






M=0 






<Z„>!t 


I 






M>0 






<Zj,>=0 




\\X 1 1 1 



FIG. 2: The bulk phase transition in the (3a — )^ plane as seen for lattice size 4 x 12'^. 

For /5f = - as one turns on the monopole suppression - the bulk transition moves down 
from (3a — 2.5 to lower /J^-values intersecting the A-axis at A ~ 0.92 as one can see from 
Fig. 121 The phases I and II are denoted as in Fig. ^ The bulk transition is characterized 
by the condensation of Z2 monopoles within phase I (M > 0) and by their suppression in 
phase II (M = 0). It has been located by monitoring the monopole density M (Fig. E)), the 
plaquette and the twist variables (for the definition of the latter see below) as a function of 
(3a for varying A. Fig. Elfor a lattice size of 4 x 12^ (i.e. for non-zero temperature) shows the 
transition to be discontinuous for small A turning into a continuous one (probably second 
order) at larger A and lower (3a, respectively (see the next paragraph). The end point of the 
first-order branch (at /5a — 1-2 for Nt = 4) seems to be identical with the lower A endpoint 
of the finite temperature transition to be reported in Section HVl For (3a > 0, A > 1.0 no 
bulk transition is observed anymore. Our /3a — A phase diagram looks very similar to that 
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FIG. 3: Monopole density as a function of (3a for different values of A (1^ = 4 x 12^) with a statistics 
of O(IO^) configurations. 



of the Villain case in the Pa — plane [22, |26[ . But note that in the Wilson case phase II 
seems to be disconnected from phase I. 

In any case the properties of the different phases should be established by means of 
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observables well defined in the continuum theory, e.g. the glueball spectrum or the screening 
mass. Such a project is currently under consideration. 



III. TWIST SECTORS AND TUNNELING 



Vortex free energies related to 't Hooft spatial loops were recently measured on the lattice 

3Q- 



within the fundamental and the adjoint representation of SU{2) |19|, |20|. In the 5*0(3) case 
the Villain discretization was used and it was stressed that on the torus T'^ the following 
equivalence holds 

t.s. crp=±l c 

where the left hand side is summed over all the twist sectors. To obtain the equivalence of 
this modified 5*0(3) Villain model with the SU (2) Wilson action case with standard periodic 
boundary conditions the global constraint 

N,,= II ap = +l. (7) 

P € plane fiu 

has to be satisfied. It is straightforward to see that an observable able to distinguish between 
trivial and non-trivial twist sectors is given by 

^t^'^ = J-J^Yl n signTrp[/p, (ep^^^ = 1). (8) 

^ per P a plane p.u 

These twist variables are again truly 50(3) observables since due to the boundary conditions 
the signs of the links in the fundamental representation drop out in the product, the plane 
extending over the whole length of the space and time directions. 

A connection between the existence of non-trivial twist sectors, the presence of center 
vortices and, in the deconfined phase, the occurence of a state characterized by a value —1/3 
for the adjoint Polyakov loop La (first observed in ^ 23]) was proposed by de Forcrand 



and Jahn j20|. The latter authors found the twist sectors to become metastable under 
local updates close to the bulk transition, i.e. the suppression of Z2 monopoles causes the 
tunnelling among different vortex vacua to become exponentially suppressed with increasing 
volume. The authors were thus prevented to go further and to investigate the possible 
extension of such connection into the finite temperature case. In the Wilson case, turning 
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FIG. 4: Monte Carlo history of the adjoint Polyakov loop and of the three electric twist variables 
at A = 0.858, /3a = 0.65 for V = 12^. 



on the 5*0(3) invariant Z2 monopole suppression term, we observed, for small volumes 
(y = 4^) and on top of the strong first order branch of the bulk phase transition (A < 0.7), 
tunneling between different distinct twist sectors 3, as well as a dynamical relation 
between the adjoint Polyakov loop and the twist observable, whose non-trivial value signals 
the presence of the state = — | also in this case. For increasing lattice volume on top of 
the bulk phase transition at low A-values the tunneling between the twist sectors becomes 
more and more suppressed. Therefore, in agreement with ref. [2^, for small values of the 
chemical potential the phase transition seems really first order as we have already argued 
earlier by monitoring other observables such as the average plaquette or the Z2 monopole 
density. 

However, as soon as the chemical potential is increased the situation changes. The ob- 
servables are no more discontinuous thus suggesting a very weak first order or second order 
phase transition. Indeed, for A > 0.75, on top of the transition, tunneling between the twist 
sectors is still evident but no suppression is seen for increasing volume as Fig. |3] shows for 
the lattice size 12^. The twists oscillate between -1 and +1 but the adjoint Polyakov l oop 
fluctuates close to zero and not between +1 and -1/3, in contrast to what was found in [20| 
for Ay = and in our case at small A. The behaviour of the twist variables reminds that of 
the fundamental Polyakov loop in the SU (2) theory in the fundamental representation, or of 
the magnetization in an Ising-like system. This justifies the definition of an order parameter 
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< z > such as 



< Z > = < -{\z^t\ + \Zyt\ + \Zzt\) > ■ 

Fig. El shows the behaviour of the susceptibihty of this new observable 



(9) 



Xtwist = iV|.((z2)_(5)2) 



(10) 



for three different volumes 8^ to 12^. A finite-size scaling test for the maximum values of 
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FIG. 5: L.h.s.: Susceptibility Xtwist as a function of A for volumes 8^,10^ and 12^ for (3a = 0.65. 
R.h.s.: Finite-size scaling of the peak values of Xtwist as a function of the linear lattice size Ns- 
The dashed line shows the finite-size scaling function fitted to the data. 



the susceptibility with the fit-ansatz j2J 



f{N,) = a{N,-\og{N,Y/'r 



(11) 



shows that the peaks increase with the lattice size with the exponent w ~ 2.4 ± .2 which 
is close to the value of the 4D Ising universality class being consistent with the theoretical 
predictions for the pure Z2 monopole system at (3a = 0. More statistics and larger volumes 
are needed to confirm this result. Moreover, the transition does not show any scaling with 
Nr, thus behaving still consistently as a T = bulk transition. 

For larger chemical potential (A > 1.0) and positive (3a, he. away from the bulk transition 
(in phase II) , tunneling is no more observed and the twist sectors become rigid. Keeping the 
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local update algorithm and turning on the Z2 monopole potential with A > 1 we can thus 
study the theory within the trivial twist sector, i.e. for the case of a zero number (modulo 
2) of (extended) vortices. 

IV. INDICATION FOR A PHASE TRANSITION AT T 7^ 

At finite temperature along a compactified dimension the periodic b.c. on the gauge 
fields allow for gauge transformations that can be classified through 7ri(S'0(3)) ~ Z2. In the 
fundamental representation this would indeed translate into a gauge transformation which 
differs by a non-trivial element of the center at the boundaries. Such transformations can 
induce a selection rule within the Hilbert space. The vacuum state, and the higher excited 
states, can exist in a superposition of the two {d + 1 for a gauge theory discretized on a 
d dimensional torus) different topological states. The creation of a vortex in the vacuum, 
as measured by the 't Hooft loop, simply means taking the expectation value between two 
vacuum states belonging to different topological sectors. A symmetry breaking argument can 
thus be applied to the vacuum state in the Hilbert space, although center symmetry is always 
unbroken. The observable z is behaving such that it averages to zero if, at finite temperature, 
the transitions occur among states which exist in a superposition of the different topologically 
distinct eigenstates, while it takes a non-zero value otherwise. 

We will identify the trivial topological sector, in which we shall be particularly interested, 
with Zxt — Zyt — Zzt — i- This corresponds to v.e.v. between the (unique even on a torus) 
topologically trivial state, i.e. no creation/annihilation of extended vortices (modulo 2). 

For what we have shown in the previous chapter, we have a local way to implement 
this constraint in the action. Indeed, since crgodicity in the — A plane is recovered 
along the bulk phase transition when it weakens, one simply needs to study the system in 
the broken phase fixing the trivial twist sector. More than a constraint it actually simply 
seems a dynamical feature of the 3-1-1 dimensional adjoint theory. Setting A > 1 and 
moving parallely to the Pa axis, we have studied the phase structure for zero {N^- — Ng) 
and non-zero temperature (A^,- <S A^^). The linear spatial lattice size was taken up to 
Ns < 16. We have determined the distribution of the fundamental Polyakov loop, governing 
the expectation value of the adjoint Polyakov loop, as an indicator for a physical temperature 
effect. Appropriate initial conditions were used to specify the trivial twist sector monitored 
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throughout the simulations. Far enough from the bulk transition the local updates are 
keeping it fixed. For sake of completeness we have studied the other twist sectors as well 
given the negative value of the Polyakov loop associated with these sectors. For A > 1 
and on the asymmetric lattice (A^,- = 4) a clear finite-T phenomenon appears. The adjoint 
Polyakov loop still averages to zero for low enough while, at higher values it eventually 
starts growing, remaining always positive (see Fig. consistently with the fixed trivial 
twist sector and with its interpretation as the exponential of the free energy of an (infinitely 
massive) adjoint quark. For increasing extent of the lattice in the imaginary time direction 
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FIG. 6: Adjoint Polyakov loop (L^) as a function of Pa in the trivial twist sector for lattice sizes 
4 X 16^ (l.h.s.) and 6 x 16'^ (r.h.s.), both for A = 1.0. 

the onset of the growth is clearly seen to be shifted to larger (3a- Alternatively one can also 
consider the spatial distribution of the fundamental Polyakov loop. For low (3a it shows a 
broad peak around Lp{x) = consistent with the Haar measure distribution ~ a/1 — L'^p, 
but with increasing (3a two symmetric peaks show up, consistent with {La) > (cf. Fig. 
The same picture occurs for Nt = 6 but with a shifted coupling as one can see from Fig. |H1 
In order to highlight this shift somewhat more quantitatively we fitted these distributions, 
produced with a statistics of 0(10^ to 10^) configurations and rescaled by the Haar measure, 
with some high order polynomials (up to 12th order). We verified that the odd derivatives are 
all zero within the errors and that the second derivative departs from zero to positive values 
(see Fig. The corresponding departure points should be taken as lower bounds for the 
critical values (3'^^^. We show them in Table I. Although we did not determine the critical 
/3a- values we see that the departure points already indicate a reasonable scaling behaviour in 
Nt. Anyway larger volumes, larger A^^-'s and also other observables will be needed to confirm 
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0.92± 0.08 


O(105) 
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1.25± 0.15 


0(10^) 



TABLE I: The /3^-values as explained in the text estimated from distributions of the fundamental 
Polyakov loop (see Fig. EJ. 

the universality with the fundamental case. Moreover, although a throughout numerical 
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FIG. 7: Spatial distribution of the fundamental Polyakov loop (Lp) in the trivial twist sector at 
A = 1.0, y = 4 X 16^ and at various /3^-values. 
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FIG. 8: Same as in Fig. d for V = 6 x 16^ 



study is still required, our data indicate that the point at which the finite temperature lines 
meet the bulk transition line coincides with the point where the bulk transition weakens 
or changes its order from 1*^* to 2"^^. This is consistent with the hypothesis that for large 
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enough volume the bulk phase transition would decouple from the deconfinement effect and 
be always weak 1*^* or 2°*^ order even at A = 0. Unfortunately, a direct numerical verification 
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FIG. 9: Second derivative of the fitting polynomial for the Haar measure rescaled fundamental 
Polyakov loop distributions at various Pa iu the trivial twist sector {V = 16^, A = 1.0); the left 
hand figure for Nf = 4, the right one for Nt = 6. 

seems unviable since from the estimates given in 0] it follows that it should occur at very 
large volumes {V ~ 0(700^) for the Villain case). 



V. CONCLUSIONS 



We have carried out a thorough study of pure SU (2) lattice gauge theory in the Wilson 
adjoint representation, decoupling the strong coupling bulk effects from the continuum limit 
by introducing a chemical potential term suppressing Z2 monopoles. As stressed in Section 
[II our formulation is not - at least not trivially - equivalent to the Villain one used in 



20|,|2ai26l in the whole parameter range we have explored. We have moreover included the 
study of the twist sectors in our analysis, which had not yet appeared in the literature when 
the previous studies in ||2^ were carried out. In this respect, we have extended the work 
of Ref. {2^ to the finite temperature case. At this point it would be interesting to study 
also the behaviour of the Villain discretization with the Z2 monopole suppression term used 
in Eq. (j2)). By monitoring mainly the twist variables and the adjoint Polyakov loop we 
studied first the character of the bulk transition. The latter turned out to become weaker 
with increasing chemical potential A turning possibly into a second order transition. At 
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strong monopole suppression A > 1 no bulk transition was found anymore. It is this area 
in the phase diagram where we started to search for a finite temperature transition. We 
have found convincing indications for such a transition at two values of A^.^ from the adjoint 
Polyakov loop and from a typical change of the distributions of the fundamental Polyakov 
loop variable. 

Due to the strong suppression of tunneling between different twist sectors the investiga- 
tions were carried out with a local update algorithm within the fixed trivial twist sector. 
Of course, a final answer should take into account also tunnelings between the sectors to be 

j30| . The results of the present investigation are collected in Fig. where we have also 
included two further points for the finite temperature transition seen for A^^ = 4 at larger 
values A = 1.5, 2.0 in an analogous way as demonstrated for A = 1.0. 
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FIG. 10: Phase diagram in the Pa — ^ plane for Nr = 4 (continuous line) and A^^^ = 6 (single 
point). 

In the meantime we are carrying out an additional investigation with the Pisa disorder 



parameter (see e.g. [3l|) adapted for the SO (3) case, which hopefully will enable us to 
determine the universality class of the transition. Preliminary data at least support the 
existence of a finite temperature transition decoupled from the bulk transition line and at 
/^A- values compatible with our results quoted here. We shall report on this in a forthcoming 
letter. 
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A quantitative study of the observed finite temperature physical transition could also 
be viable relying on pure thermodynamic quantities, as already done in I22II . through the 
Lee- Yang zeros technique or by studying the free energy of a Z2 vortex jiol . Isj ]. More careful 
investigations of this kind will hopefully tell us something about the real nature of the finite 
temperature transition we have reported here. But the mere occurence of such a transition 
is a prerequisite to a still missing proof of universality between SU{2) and 5*0(3) lattice 
gauge theories. 
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